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Stability Optimization of Laminated Composite Plates 


By 

Yoichi Hirano 


Summary. This paper is concerned with the opiimi'm design of plates with orthotropic 
layers under axial compression and shear, The plates considered are the laminates of N 
ortholropic layers whose principal material axes coincide with the plate axes, Each layer is 
assumed to have the same thickness and an equal number of fibers in the direction of +ot, 
and —a., with respect to the plate axis, The fiber directions which give the highest axial 
buckling stress and the highest shear buckling stress arc found by utilizing a mathematical 
optimization technique for various aspect ratios of the plates. Inhomogeneity in the direction 
of the plate thickness (stacking sequence) is taken into account in this analysis, 
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= plate length in the xr-direction (Fig. 1) 

= extensional stiffness of a laminated plate 
= plate length in the >’-direction (Fig. 1) 

= coupling stiffness of a laminated plate 
= bending stiffness of a laminated plate 

= Young’s moduli of a unidirectional composite parallel and 
transverse to the directions of the fibers, respectively 
= shear modulus of a unidirectional composite 
= thickness of each layer 
= Eqs. (17) 

= the ratio of N y to N x 

= number of half waves in the .v- and j’-directions, respectively 
= number of layers 

= applied force per unit length in the x- and ^-directions, respectively 
= applied shear force per unit length 
= reduced stiffness 
= transformed reduced stiffness 
= Eqs, (14) 

= Eqs. (8) 

= displacements in the a--, y- and --directions, respectively 
= displacement amplitudes defined by Eqs. (10) 

= Eq. (18) 

== absolute value of fiber di. actions with respect to x in the ;-th layer 
=; Eqs. (12) 

= Poisson’s ratio of a unidirectional composite 
=s nondimensional critical buckling stress defined by Eq. (9) 
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inverse o!' nondimensional critical buckling stress; defined by Hq, 
(16) 


I, INTRODUCTION 

Recently, lilainentary composite materials have been suggested for primary 
structures of aircrafts and spacecrafts. The reason is mainly due to the weight savings 
which can be attained. There are many design criteria in applying the composite 
materials to structures. One of them is the buckling criterion. Several theoretical and 
experimental papers 1 ’ have been published on (tie buckling of laminated composite 
plates under axial compression and shear. Most of them give their results only for such 
special cases as angle-ply or cross-ply plates, Therefore, we do not have enough 
information to design the laminated composite plates. 

This paper will present a method to design the laminated plates (Fig. 1) with 
orthotronic layers under uniaxial or biaxial compression and shear. The design 
criterion is the buckling stress. Each layer of the plate is assumed to have the same 
thickness and an equal number of the same kind of libers in the + a, and —a, directions 
with respect to the .v coordinate in the same type of matrix. Therefore, each layer can be 
considered to be orthotropic. Inhomogeneity in the direction of the thickness of the 
plate (stacking sequence) is taken into account in the calculation, 

The present problem is to find the fiber directions of all the layers that give the 
highest buckling stress and, therefore, is an unconstrained maximization problem. The 
objective function and the design variables are the critical buckling stress and the liber 
directions respectively. Preassigned parameters are the material properties, the 
thickness of each layer, the number of layers, and the aspect ratio of the plates. The 
optimization technique used is Powell’s method (conjugate direction technique). This 
method is one of the best ones to find the optimum without using the derivatives of the 
objective function. 

2. Derivation of Buckling Differential Equations 

Extensional, coupling and bending stillness, which are expressed as and Dy 

respectively, are first introduced. They are defined in terms of transformed reduced 


stillness as follows. 

(zl U . By, i)y) « f (? 0 (1. = , = 2 * )<h (I) 

J 1.2 

The transformed reduced stillnesses of each layer are calculated by the following 
equations. 

0it vos 4 oCj + 2 (Q l2 -I- 2()„ {1 ) siira, cos 2 a, + Q 21 sin 4 a, (2a) 

Q\z = (0,t + Qzz ~ 4(?(,„)sin-a,.cos 2 !X i + Q l2 (sin 4 a f + cos 4 a,) (2b) 

Qn “ 0usin + a, + 2 (Q l2 + 2() t , (> ) siira, cos 2 ot, + Q 22 cos 4 « f (?o) 

0th = 0 (2d) 
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Qu, - o 

Q m = «?1 1 + Qu ~~Qis " 2(?„ t ,)sin 2 «,cos J a, + <?<, ( ,lsin 4 a, + cos A x,) 
where Q t j s are reduced stillnesses and tire defined as 

<?n = K|.'(l -Vuv,,) 

Q 1 1 r [ i t ( 1 ■“ 

= 15, (I - v, 2 v 2l ) 


(2e) 

(20 


(3a) 

(3bl 

(3c 


/) - (.HI 

Vd(. — ci|.j v 

In Bqs. (2) (5 lt) and (j 2b are equal to zero, because the number of libers in the -I a, and 
-a, direction, are the same, lixtensional, coupling and bending stillnesses are 
calculated for the present ease as 
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where /, is ills' thick™ or *1. layer. N is «« *»*“• ol ' “ ra ' 

subscript of is thc nuu ' ber 0 | ° aC '3;. cmnti o„s for the general laminated 

i , , t i A — /}. , U*, vv \ ' (Hi: 'h 
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,i„ r ,he present ease ft. in Be,. ,M is et,ua, .0 «r„. a,ul the M*. ^.bnnatious 
arc assumed as 


i, |7 cos (jun.x. it) siu (imyh) 
r=r sin (nmxiti) cos (nny;b) 
w=-.\y sin {nmxiti) sin (»ny b) 


(6a) 

(6b) 

(6c) 


These ^relations satisfy the simply 

,,^ 0i /, ( S2 of Ref. 9). Substituting bqs. (6) mto l.qs. (M ano t > 
following buckling formula. 
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T\ 2 = (A n + A bb )tnim 2 (a/h) (8b) 

r l3 = By ,mV + (By 2 + 2B bb )intrn 3 (a/b) 2 (8c) 

T '12 — A bb m 2 n 2 + A 2 1 n 2 n 2 (a/b) 2 (8d) 

T' Z3 = (By 2 + 2B bb )m 2 nn 3 (a/h)+B 22 n 3 n 3 ((i/h) 3 (8e) 

r 3 , = /3, ,ni 4 tt' i + 2 (D 1 2 + 2D bb )m 2 n 2 n\a/h) 2 +• D 22 n 4 n%ifb) 4 (8f) 

To get the critical hackling stress the smallest value of Eq. (7) must be found by a 
searching procedure involving integer values of rn and n. The critical buckling stress is 
denoted by (NJJt and a new notation </> is introduced, 

</' = l2(^JJr/(n 2 t 3 Q 22 ) ( 9 ) 

For isotropic plates (ft is equal to 4, when alb — I and A'=0, 


4. Calculation oh Shear Buckling Stress 

For the present case N x and N y are equal to zero, and the plates are assumed to be 
simply supported (S2 of Almroth) at four edges. The following deflection function 
satisfy the boundary conditions. 


or • 


v-, ^ _ nmx , tmv 
" = I L u m„cos — - • sin 

III :■ I II =1 a " 

_ . nmx tmv 

l '= L L i, sin — cos 

in = I ii = i 11 ” 

_ . nmx , tmv 

»'= L L »' sin -- - sin , 

III =: I II = 1 « O 


Substitution of Eqs. (10) into Eqs. (5a) and (5b) yields: 

««ii«„„, + ««i2‘ 7 „,„ “ tta 13 w,„ n = 0 

«*12'7,i,„ + act 21 f„ w - na 2i w nl „ = 0 


where 


a i i — Ayyttr + A bb n~R 2 
= (^i 2 + A bb )mnR 
a, 3 = By yin 3 + (B 12 + 2B bb )mii 2 R 
ot 22 = A bb tn 2 + A 22 trR 2 
a 22 ~ (B l2 + 2B bb )tn 2 ttR + B 22 ii 3 R : 


(10a) 

(10b) 

(10c) 

(11a) 

(lib) 

(12a) 

(12b) 

(12c) 

(12d) 

(12c) 


In the above expressions R is the aspect ratio alb of the plates. From Eqs, (11a) and 
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(lib) tl nm and tan be expressed in terms of 
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~ *13*22 — *12*23 

(14a) 

7(iiii — *11*23 ~ *12*13 

(14b) 

— «i 1*22 ~~ *1 2~ 

(14c) 


Now Eqs, (13) arc substituted into Eq, (5c) and Galerkin's method is applied. Then we 
get the following expression. 

t t II * = 0 (15) 

* pq m e= 1 II - 1 

where /; and q arc the number of half waves in the x- and ^'-directions respectively; V', 
H m „ and V m arc defined as: 


W xr 32irR 

II„,„ — --?■ ~ > i ~ ~~i when /> + m: odd and q + n: odd 
p- — nr q~ — ir 

= 0 when p + nt: even or q + n: even 

= {»n/> 4 + 2(Z)|, + 2 D ( Jp 2 q 2 R 2 + 0w/*K 4 } 

- {BnP 3 + ( 0.2 + 2B M MfR 2 }S„,/U lul 

- {(B |2 + 2 B ( Jp 2 qR + B 21 q>R>}Cr tul /U, ul ) 


(16) 


(17a) 

(17b) 


(18) 


Eq. (15) is a system of homogeneous linear equations in ir,„„. This system can be divided 
into two groups, one containing for which nt + it are odd and the other for which 
»/+« are even. Two buckling forces are obtained from these two groups and the lower 
one corresponds to the critical buckling force of the laminated plates. 

A computer program to solve Eq. ( 1 5) was written, and was checked for the case of 
isotropic plates by comparing with the results obtained by Stein and Neff [ 70], 


5. Method of Optimization 


The problem is to find the fiber directions which give the maximum critical buckling 
stress without any constraints. Therefore we can apply one of the unconstrained 
optimization techniques. Since the objective function for the case of axial compression 
is the rather complicated function of the design variables and the objective function for 
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Table 1, Material properties 




Property 


Boron/Epoxy 





Ei 

*'|2 
O’, 2 


2.11 x It) 4 kg/mm 2 (30 x 10° psi) 
2.11 x 10 3 kg/mm 2 (3x 10 6 psi) 
0.3 

7.03 x I0 2 kg/mm 2 (1 x 10 fl psi) 



Table 2. Optimum fiber directions for 3-laycrcd plates 
(<j,b~- 1, k-0) under axial compression 


Cusc 



Fiber directions 
(in degree) 

Critical 

stress 

Number of 
half waves 

«i 

a 2 


<l> 

m n 


s 

0.0 

0.0 

0.0 

12,921 

i i 

i 

F 

45.0 

45,0 

45.0 

22,000 

i t 


s 

0,0 

90.0 

0.0 

12,921 

i i 


P 

45.0 

135.0 

45.0 

22,000 

i i 


s 

45.0 

0.0 

45.0 

21.664 

i i 


F 

45.0 

0.0 

45.0 

21.664 

i i 


S 

0.0 

45.0 

0.0 

13.25S 

i i 


F 

45.0 

44.9 

45.0 

22.000 

i i 


S 

30.0 

30.0 

30.0 

19.730 

i i 


F 

45.0 

45.2 

45,0 

22.000 

i i 


S 

90.0 

0.0 

90.0 

9.671 

2 1 

6 

F 

135.0 

45.0 

135.0 

22.000 

1 1 


S 

45.0 

45.0 

45.0 

22.000 

1 1 

7 

F 

45.0 

45,0 

45.0 

22.000 

1 1 


S: Starting values, F: Final optimum values. 


the case of shear is not obtained explicitly, optimization methods without using the 
derivatives are suitable for solving the problem, Powell’s method [7 7] (conjugate 
direction technique) is selected for use, since it is one of the best methods to find the 
optimum without using the derivatives [72], This method is intuitively explained by 
Fox [75] as follows: “Given that the function has been minimized once in each of the 
coordinate directions and then in the associated pattern direction, discard one of the 
coordinate directions in favor of the pattern direction for inclusion in the next 
minimizations, since this is likely to be a better direction than the discarded coordinate 
direction. After the next cycle of minimizations, generate a new pattern direction and 
again replace one of the coordinate directions,” 

Powell's method is now applied to find the maximum value of tj) and the minimum 
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Tabu; 3. Optimum liber directions lor 4-layered plates 
(«//>=!, k — 0) under axial compression 


Case 



Fiber directions 
(in degree) 


Critical 

stress 

<l> 

Number of 
half waves 

#i 

»2 


«4 

m n 


S 

0.0 

0.0 

0.0 

0,0 

12.921 

i i 


F 

45.0 

45.1 

45.0 

45.0 

22,000 

i i 


S 

30.0 

30.0 

30,0 

30.0 

19.730 

i i 


F 

45.0 

45.1 

45.1 

45.0 

22.000 

i i 


S 

45.0 

45.0 

45.0 

45.0 

22.000 

i i 


F 

45.0 

45.0 

45.0 

45,0 

22,000 

i i 


S 

90,0 

90.0 

90.0 

90.0 

8.421 

2 1 


F 

135.0 

135.0 

135.0 

135.0 

22.000 

i i 


S 

90.0 

0.0 

0.0 

90.0 

12.640 

2 1 


F 

135.0 

45.0 

45.0 

135.0 

22,000 

1 1 


S 

45.0 

0.0 

0,0 

45,0 

20.865 

> 1 


F 

45.0 

45.1 

45.0 

45.0 

22.000 

1 1 


S 

0.0 

45.0 

45.0 

0.0 

14.056 

1 1 


F 

45.0 

45.0 

44.9 

45,0 

22.000 

1 1 


S 

10.0 

20.0 

30,0 

40.0 

14.431 

1 1 


F 

45.0 

45.0 

44.9 

45,0 

22,000 

1 1 


S: Starting values, F: Final optimum values, 


value of Starting values of fiber directions (a,, a 2 • > • , a„) are necessary to begin the 
calculation, and the new fiber directions which give the higher buckling stress are 
obtained after each iteration, Powell’s method requires that the objective function be 
unimodal. But we do not know if the function is unimodal or not. Therefore, trials with 
several starting points are desirable, 

6. Numerical Results for the Case of Axial Buckling 

Numerical calculations were made for the laminated plates with three, four and six 
layers. The plates considered have the various aspect ratios and are under uniaxial or 
biaxial compression. Seven or eight different combinations of fiber directions are used 
to start the calculation. Computer code developed by Powell was combined to the code 
written for the present problem. 

Convergence limits for all the design variables were set to be equal to 0.1 L and 
maximum step size multiplier [14\ in single variable search was set to be equal to 10.0. 
The materials considered are Boron/Epoxy and the properties [5] are shown in Table 1 . 
The thickness of each layer is assumed to be 0.254 mm (0.01 in.). The buckling formula 
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Tahi.h 4. Optimum liber directions for 6-layered plates 
(«//>=!, A*= 0) under axial compression 


Case 

— 



Fiber directions 
(in degree) 



Critical 

stress 

'l> 

Number of 
half waves 



“j 

«4 

«5 

<*:. 

/u n 


s 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

12.921 

i i 

1 

F 

45.0 

45.1 

45.1 

44.9 

45,1 

45.0 

22,000 



S 

30.0 

30,0 

30.0 

30.0 

30,0 

30.0 

19.730 

i i 


F 

45,0 

45.0 

45.0 

44.9 

45,0 

45.0 

22.000 

i i 


S 

45,0 

45.0 

45.0 

45,0 

45.0 

45.0 

22.000 

i i 

.5 

F 

45.0 

45.0 

45.0 

45.0 

45.0 

45.0 

22.000 

i i 


S 

90.0 

90.0 

90.0 

90.0 

90.0 

90,0 

8.421 

2 1 

4 

F 

135,0 

135.0 

135,0 

135.1 

134,9 

135.0 

22.000 

1 1 


S 

90.0 

0.« 

90 0 

0,0 

90,0 

0,0 

12.272 

1 1 

j 

F 

135.0 

4 f j! 

135.4 

44.8 

135.0 

45,0 

22,000 

1 1 


S 

45.0 

4,'j.O 

0.0 

0.0 

45.0 

45.0 

21.664 

1 1 

6 

F 

45.0 

45.0 

45.0 

45.1 

44,9 

45.0 

22,000 

1 I 


S 

0.0 

0.0 

45,0 

45.0 

0.0 

0.0 

13,258 

1 1 

7 

F 

44.9 

44.9 

47.4 

45.4 

45.1 

45.1 

21.999 

1 1 


S 

10.0 

20.0 

30.0 

40,0 

50.0 

60,0 

11.711 

1 1 

0 

F 

45.0 

45.0 

44.8 

45.3 

44.6 

44.9 

22,000 

1 1 


S: Starling values, F: Final optimum values, 


is a function of half waves in the x- and j’-directions. Therefore, the numbers of half 
waves in the x- and ^-directions were varied from 1 to 10 and from 1 to 5, respectively, 
to get the buckling stress for the assigned fiber directions, 

The results for three- and four-layered plates with ajb— 1 and A'=0 are presented in 
Table 2 and Table 3, respectively. These two tables show that the results obtained do 
not depend on the starting values of fiber directions except in Case 4 of Table 2. This 
case shows that the calculation converged to local maximum, To show the numerical 
convergence two examples are shown in Figs. 2 and 3 for Case 4 of three-layered plates 
and for Case 8 of four-layered plates, respectively. In these figures the abscissa is the 
number of iteration, Each iteration includes many function evaluations. From these 
figures and tables it can be concluded that this method of finding the best fiber 
directions works well. Then, the method was applied to six-layered plates and some of 
the obtained results are shown in Tables 4, 5 and 6. Table 4 is for the case of ajb= 1 and 
k = 0. This table shows that the present method also works well for six-layered plates. 
Almost all the fiber directions obtained are close to 45°, but some of the directions are 
not close to 45 because of the slow convergence of the numerical calculations. Table 5 
is for the case of a/b= 0.5 and 0.5. The table shows that the final critical buckling 
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Taiu r 5, Optimum fiber directions for 6-layered plates 
(</,/> e- 0,5. A- -0.5) under axial compression 


Case 




Fiber directions 
(in degree) 



Critical 

stress 

1> 

Number of 
half waves 



*3 

*4 

«5 


w n 


s 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

34.457 

1 2 


F 

13.2 

0.0 

0.0 

0.2 

15.9 

0,0 

37.011 

1 1 


S 

30.0 

30,0 

30. ,i 

30.0 

30.0 

30,0 

32,633 

1 1 


F 

3.0 

11.7 

35.7 

-21.9 

0.1 

11.6 

36.966 

1 1 


S 

45.0 

45.0 

45.0 

45.0 

45,0 

45.0 

26.016 

1 1 


F 

4.0 

10.0 

36.1 

-2.0 

-21.8 

2.1 

36.934 

1 1 


S 

90.0 

90.0 

90.0 

90,0 

90.0 

90.0 

7.486 

1 1 

4 

F 

186.3 

179.3 

90.5 

180.0 

189,3 

167.9 

36.492 



S 

90.0 

0.0 

90.0 

0.0 

90.0 

0.0 

21.294 

1 1 


F 

173.2 

0.2 

146.6 

0.1 

157.2 

1.8 

36.919 

1 1 


S 

45.0 

45.0 

0.0 

0.0 

45.0 

45.0 

26.441 

1 1 


F 

0.0 

2.4 

138.7 

0.4 

-25.3 

2.4 

36.796 

1 1 


S 

0.0 

0.0 

45.0 

45.0 

0.0 

0.0 

35.353 

I 1 


F 

10.7 

-0.9 

36,4 

39.1 

76,9 

4.9 

36.909 

1 1 


S 

10.0 

20.0 

30.0 

40.0 

50,0 

60,0 

18.566 

1 1 

C> 

F 

-7.9 

-3.8 

22, 1 

-0.5 

-7.4 

12.3 

37.024 

1 1 


S: Starting values, F: Final optimum values, 


stresses obtained for the different starting values are almost the same but the 
corresponding fiber directions are not the same. This may be due to the fact that the 
objective function is not unimodal for this case. Table 6 is for the case of a/b= 1.0 and 
A = 0.5. Summary of the numerical results is given in Table 7. In this table the rows with 
an asterisk show that the final results obtained depend on the starting values and the 
values shown correspond to the highest critical buckling stress obtained among 8 cases. 
In these cases the critical buckling stresses obtained are not much different from each 
other, but the fiber directions highly depnd on the starting values. The fiber directions 
in the rows without an asterisk have no decimal, because almost all the directions 
obtained for seven or eight cases are close to the values shown. 

The computer used was IBM 360/67 and average epu time to calculate eight cases for 
a six-layered plate under k = 0 was 158 seconds. 
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Table 6, Optimum fiber directions for 6-laycred plates 
Ui/h— I, A- *0.5) under axial compression 


Case 




Fiber directions 
(in degree) 



Critical 

stress 

<I> 

Number of 
half waves 

«i 



« 4 


»6 

m n 


s 

El 


mm 


1 


8,614 

i i 


F 




44,8 


45.0 

14,667 

i t 


S 

30.0 

30.0 

30.0 

30,0 

30.0 

30.0 

13.154 

i i 


F 

45.0 

45.0 

45,3 

45,1 

44.9 

45.0 

14.667 

i i 


S 

45.0 

45.0 

45.0 

45.0 

45.0 

45,0 

14,667 

i i 


F 

45.0 

45.0 

45.0 

45.0 

45.0 

45.0 

14.667 

t i 


S 

90.0 

90.0 

90.0 

90,0 

90.0 

90.0 

7.486 

2 1 


F 

135,1 

134,9 

134.9 

135.0 

135.0 

135.0 

14.667 

1 1 


S 

90,0 

0,0 

90,0 

0,0 

90,0 

0,0 

8.182 

1 1 


F 

134.9 

44.8 

134.9 

44.9 

134.9 

45.0 

14,667 

1 1 


S 

45.0 

45.0 

0.0 

0.0 

45.0 

45.0 

14.442 

1 1 


F 

45,0 

45.1 

44.7 

44,9 

45.1 

45.0 

14.667 

1 1 

7 

S 

0,0 

0,0 

45.0 

45.0 

0.0 

0.0 

8.838 

1 1 


F 

45.0 

45.0 

45.1 

45,1 

45.1 

45,0 

14.667 

1 1 

Q 

S 

10,0 

20.0 

30.0 

40.0 

50.0 

60,0 

7.807 

1 1 


F 

45.0 

45.1 

44.7 

45.3 

45.0 

45.0 

14.667 

1 1 


S: Starting values, F: Final optimum values. 


7. Numerical Results for the Case of Shear Buckling 

Powell’s computer code was rearranged into a code for the calculation of the shear 
buckling stress. The material considered is Boron/Epoxy (Table 1), The thickness of 
each layer is assumed to be 0,254 mm (0.0 1 in). The buckling stresses were calculated by 
taking w=l~3, «=1~3 for /?=1; and /h=1~5, h= 1 ~ 5 for R=l,5 and 3. The 
numerical errors of calculated buckling stresses for R= l and 1.5 are less than 3% and 
the error for R—3 is 1 1%, when all fibers are in the direction of 90 1 ’ with respect to the 
,v-axis. Therefore, the obtained results for R= 1 and 1.5 are accurate enough, but the 
results for /?= 3 may not be accurate enough. 

Numerical calculations were first made for the case of three-layered plates with /?= 1 
and an example of the numerical convergence is shown in Figure 4. In this figure the 
abscissa is the number of iterations, and each iteration includes many function 
evaluations. Then calculations were made for the case of six-layered plates with R = 1, 
1.5 and 3. The results for these cases are presented in Tables 8, 9 and 10. Asterisks in 
these tables indicate that the numerical calculation was stopped because of the fact that 
a maximum change in a single variable search did not alter the objective function value. 
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Table 7, Summary of optimum fiber directions for the cuse of axial compression 


No , or 
layers 

k 

alb 

«l 

«2 

«3 

«4 

“5 


<l> 

3 

0 

1.0 

45 

45 

45 " 




22,000 

4 

0 

1,0 

45 ' 

45 

45 ' 

45 ‘ 



22,000 

6 

0 

0.5 

0 

0 

0 

0 

0 

0 

42.171 

6 

0 

0,8 

38 ' 

38 

38 

38 ' 

38 ' 

oo 

23.154 

6 

0 

1.0 

45 

45 ' 

45 ' 

45 

45 

45 * 

22,000 

*6 

0 

1,25 

49 . 9 ' 

51.0 

48,6 

48 .rf 

51.0 

49.9 

23.116 

6 

0 

2,0 

45 

45 

45 

45 ' 

45 

45 

22.000 

*6 

0,5 

0,5 

7,9 

2.8 

22.1 

0.5 

7.4 

12.3 

37.024 

6 

0.5 

1.0 

45 

45 

45 

45 

45 ' 

45 

14.667 

♦6 

0.5 

2.0 

67.1 

56.4 

56,2 

55 . 5 ' 

64.0 

61.4 

12.556 

6 

1.0 

1.0 

45 ' 

■■ 

■■ 



45 


♦6 

1.0 

2,0 

71.6 

HH 

E9 



74. i 1 

8.051 



No. of iterations 


Fig. 4. Variation of i, and i p with number of iterations. 
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Tam i 8. Optimum fiber directions for O-layered pltilcs (c/,-7 1 I) muler shear 


‘use 




l ; iber directions (in degree) 



\ uv 




v *1 


■X. 





s 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0753 

1 

i-' 

45.6 

0.0 

0.2 

0.1 

0,0 

43.7 

1 .2708 


s 

30.0 

30.0 

30,0 

30.0 

30.0 

30.0 

1.2665 

** 

I- 

45.2 

44.0 

45.0 

45.1 

45.1 

45,1 

1.1283 

3 

s 

45.0 

45.0 

45.0 

45.0 

45.0 

45.0 

1.1283 

1’ 

45,0 

45.0 

45.0 

45.0 

45.0 

45.0 

1.1283 

4 

s 

00.0 

0.0 

00.0 

0.0 

00.0 

0.0 

1.00 II 

1’ 

135.0 

45,0 

134,8 

44.4 

134.0 

45.0 

1.1283 

5 

s 

F 

45.0 

45.0 

0.0 

♦ 

0.0 

45.0 

45.0 

1.1431 


s 

0.0 

0.0 

45,0 

45.0 

0.0 

0.0 

1.0103 

(> 

F 

44.0 

45,3 

45.4 

45.3 

45.2 

45.0 

1.1283 


S 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

2.0520 

7 

V 

45 1 

45.1 

44,3 

45.1 

45.0 

45.0 

1.1283 

8 

s 

1' 

0.0 

0.0 

4M) 

♦ 

45.0 

00.0 

00.0 

3 4567 

a 

s 

I' 

00.0 

00.0 

oo.o 

* 

00.0 

00.0 

00.0 

1.0753 

10 

s 

60.0 

60.0 

60.0 

60.0 

60.0 

60,0 

1.2665 


45.1 

45.0 

43.8 

45.0 

45.1 

44.0 

1.1283 


S: Starting values, !■': I’inal optimum values. 


The computer used for the present case was IBM 3033 and eptt time to obtain Table 0 
wtis 1087 seconds. 


8. Conclusions 

A method to linn the best liber directions of the laminated plates under axial 
compression and shear lias been presented in this paper. From Table 7 it can be said 
that the best fiber angles in all layers are 45 for the plate with «//>-- 1 and 2 under 
uniaxial compression. For the plates with «;/>•-=(), 5, 0.8 and 1.25 under k 0 the best 
angles in all layers are 0‘ . 38' and 50 respectively. It is interesting to note that the best 
fiber directions lor the case of k~ 0 is the same in all layers, even if the stacking 
sequence is taken into account. For the plates under A ~ 0, simple conclusions cannot 
be obtained. From Tables 8, 9 and 10 the following conclusion can be obtained for the 
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Tabi.p. 9. Optimum fiber directions for 6-laycred plates (ajb= 1.5) under shear 


Fiber directions (in degree) 

Case \<i x I0 5 





«2 


«4 

“s 



I 

S 

0.0 

0,0 

0,0 

0.0 

0.0 

0.0 

12.279 

1 

F 

54.6 

54.5 

54.4 

54.3 

54,6 

54.5 

4.4804 


S 

30,0 

30.0 

30.0 

30.0 

30,0 

30,0 

6.0322 

z 

F 

54.7 

54.4 

54.6 

54.6 

54.7 

54,5 

4.4805 

T 

S 

45,0 

45.0 

45.0 

45.0 

45,0 

45.0 

4.6801 


F 

54,5 

54.7 

54,2 

54.3 

54.8 

54.4 

4,4805 

4 

S 

F 

90.0 

0,0 

90,0 

* 

0,0 

90.0 

0.0 

6.8648 

5 

S 

F 

45.0 

45,0 

0.0 

* 

0.0 

45.0 

45.0 

4.7749 


S 

0,0 

0.0 

45.0 

45.0 

0,0 

0.0 

11.516 

0 

F 

54.5 

54.5 

55.2 

55.1 

54.4 

54.7 

4.4805 

n 

S 

10.0 

20.0 

30.0 

40,0 

50.0 

60,0 

8 3499 

/ 

F 

54.6 

54.3 

54.6 

54.1 

54.6 

54.5 

4.4805 

O 

S 

0,0 

0.0 

45.0 

45.0 

90.0 

90.0 

12.486 

O 

F 

54.4 

54.7 

54.1 

54.6 

125.6 

125.5 

4.4805 

n 

S 

90.0 

90.0 

90.0 

90.0 

90.0 

90.0 

6.7886 

V 

F 

126.0 

125.4 

90.3 

90.4 

124.9 

126.3 

4,5004 

in 

S 

60.0 

60.0 

60,0 

60.0 

60.0 

60,0 

4.5425 

1 u 

F 

54.6 

54.4 

54.3 

54.1 

54.4 

54.6 

4.4805 

S: 

Starting vi 

dues, F: 

Final optimum 

values. 






case of shear buckling. An angle-ply laminate gives the highest shear buckling stress, 
even if almost complete freedom is given in the selection of fiber directions. The best 
fiber directions for the cases R= 1, 1.5 and 3 are 45\ 55" and 60" respectively. These 
angles are equal to the ones obtained by Housner-Stein for the case of angle-ply 
laminated plates. 

This work was done during the author’s stay at Rensselaer Polytechnic Institute as a 
visiting associate professor and was supported by NASA under Grant No. NGL-33- 
018-003. The author wishes to acknolwcdge the helpful advice of Professor Nicholas J. 
Hoff. 
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Tabi i 10. Opliimim liber directions for 6-lnyered plates (a'h- 3) under shear 


Fiber directions (in degree) 

Case -- ft, x | ()" 





«a 

“4 


«o 


s 

F 

0.0 

0.0 

0.0 

* 

0.0 

0.0 

0,0 

23.473 

S 

30.0 

30,0 

30,0 

30.0 

30.0 

30.0 

11.261 

F 

60.6 

61,4 

61.5 

66.9 

60,8 

59.6 

6.5480 

S 

45,0 

45.0 

45.0 

45.0 

45.0 

45.0 

7.5115 

F 

61.2 

58.5 

62.3 

62.7 

61.0 

59.6 

6.5484 

s 

90.0 

0.0 

90,0 

0.0 

90.0 

0,0 

10.015 

F 

90,0 

35.4 

90.0 

0.0 

90.0 

0.0 

9.6774 

S 

45.0 

45.0 

0.0 

0.0 

45,0 

45.0 

7,6926 

F 

60.6 

60,6 

58.5 

57,9 

60.7 

60.5 

6,5469 

S 

0.0 

0.0 

45.0 

45.0 

0.0 

0,0 

21,521 

F 

60,7 

60.5 

62.3 

66,9 

59.8 

60.2 

6.5477 

S 

10.0 

20.0 

30,0 

40.0 

50,0 

60.0 

13,175 

F 

60.9 

60.8 

60,1 

69.2 

59,8 

60.0 

6.5484 

S 

0.0 

0.0 

45.0 

45.0 

90,0 

90.0 

18.204 

F 

59.3 

121.0 

88,2 

58.4 

116,8 

1 19.5 

6.5558 

S 

90.0 

90.0 

90.0 

90.0 

90.0 

90.0 

8.2017 

F 

90.0 

90.0 

90.0 

90,0 

90.0 

90.0 

8.2017 

S 

60.0 

60.0 

60.0 

60.0 

60.0 

60.0 

6.5478 

F 

60.7 

60.0 

60.0 

60.0 

60.8 

60.4 

6.5476 


S: Starting values, Ft Final optimum values. 
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